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ABSTRACT 


The  development  of  ne%f  decision  support  systems  for 
Antisubmarine  warfare  will  entail  the  installation  of 
propagation  loss  models  on  ASW  aircraft.  The  decision  to  put 
either  a  range  dependent  or  range  independent  model  in  the 
system  will  affect  the  predicted  ranges,  the  overall 
probability  of  detection,  and  the  computation  time. 

Comparisons  of  the  range  dependent  ASTRAL  and  range 
independent  RAYMODE  propagation  loss  models  were  made  in  the 
Eastern  Mediterranean,  the  Gulf  of  Oman  and  the  South  China 
Sea  for  eight  source/receiver/freguency  combinations. 
Computation  time  differences  between  the  two  models  were  not 
significant  at  either  of  the  source  frequencies  (50  Hz  or  400 
Hz).  RAYMODE  Showed  much  better  correlation  with  the  split 
step  PE  model  which  was  used  as  a  standard.  The  ASTRAL  model 
often  predicted  lower  transmission  losses  than  either  RAYMODE 
or  PE.  For  the  short  detection  ranges  normally  encountered  in 
air  ASW  the  more  complex  range  dependent  models  are  not 
necessary.  The  RAYMODE  model  or  a 
independent  model  will  provide  adequate 
predictions . 
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1 .  INTRODUCTION 


The  intent  of  this  paper  is  to  present  a  comparison  study 
of  the  effect  of  using  a  range  dependent  versus  a  range 
independent  propagation  loss  model  in  an  airborne  ASW 
tactical  decision  support  system.  In  an  ASW  tactical 
situation,  processing  time  and  accurate  probability  of 
detection  are  extremely  important.  The  ability  to  produce 
guick  and  reliable  transmission  loss  information  while  on 
station  would  be  a  definite  asset  in  a  decision  support  role. 
The  two  transmission  loss  models  currently  being  considered  as 
candidates  to  be  installed  in  an  aircraft  decision  support 
system  are  the  range  dependent  ASTRAL  or  the  range  independent 
RAYMODE  model  [Ref.  1].  Using  a  range  dependent  (RD)  versus 
a  range  independent  (RI)  propagation  loss  model  may  affect 
both  the  processing  time  and  the  probability  of  detection. 

This  study  compares  the  ASTRAL  and  RAYMODE  models  in 
several  areas.  First,  it  will  review  the  physics  that  both 
models  use  to  develop  their  propagation  loss  information. 
Next,  using  three  different  geographical  areas  and  a 
combination  of  two  source  frequencies,  two  source  depths,  and 
two  hydrophone  depths,  transmission  loss  profiles  will  be 
generated.  The  profiles  are  then  used  to  develop  detection 
ranges  and  probabilities  of  detection  for  three  specific 
sonobuoy  patterns  (2x8,  4x4,  5x6x5,  shovm  in  Fig.  1,  spacing 
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is  based  on  the  predicted  detection  ranges) .  The  conparisons 
of  the  produced  transmission  loss  profiles,  the  detection 
ranges,  and  the  probabilities  of  detection  of  the  patterns 
will  all  be  presented.  The  processing  time  required  by  each 
model  to  run  the  transmission  loss  profiles  will  also  be 
discussed . 


Figure  1  -  Sonobuoy  Pattern  Geometries 
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IZ.  PROPAGATION  LOSS  MODELS 


The  Navy  currently  uses  many  different  oceanographic  and 
propagation  loss  models.  Propagation  loss  is  the  attenuation 
of  the  sound  intensity  as  it  passes  through  the  ocean  from  the 
source  to  the  receiver.  Propagation  losses  are,  in  general, 
affected  by  such  factors  as  the  temperature /depth  profile  of 
the  water  column  from  the  surface  to  the  bottom,  the  source 
frequency,  the  condition  of  the  ocean  surface  (wind  speed  or 
sea  state) ,  bottom  depth,  bottom  composition  and  thickness, 
salinity,  source  and  receiver  depths,  and  range  to  the  source. 
Once  a  model  is  tested  and  accepted  by  the  Navy  it  is 
designated  as  a  Navy  Standard  Model.  There  are  both  range 
dependent  and  range  independent  standard  models  .  The  range 
independent  models  assume  the  water  column  does  not  change 
with  range  and  do  not  take  into  account  the  ocean  bottom 
topography.  A  single  sound  speed  profile,  water  depth  and 
bottom  loss  are  used  for  the  calculations  for  an  entire  area. 
The  range  dependent  models  do  use  bottom  topography  and  allow 
the  sound  speed  profile  and  bottom  loss  to  change  with  range. 
RAYMODE  is  currently  the  only  range  independent  Navy  Standard 
model,  while  ASTRAL  and  PE  (Parabolic  Equation  model)  are  the 
most  commonly  used  RD  models.  The  RAYMODE  and  ASTRAL  models 
use  distinctly  different  methods  for  calculating  the 
transmission  loss  from  a  source  to  a  receiver.  In  general. 
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most  propagation  loss  models  rely  on  either  wave  theory  or  ray 
theory  [Ref.  2:p.  1]. 

The  wave  equation  is  a  partial  differential  equation 
relating  the  acoustic  pressure  p  to  the  coordinates  x,y,z  and 
the  time  t,  and  may  be  written  as 


dt^  dx^  dy^  dz^ 


(1) 


where  c  is  the  sound  speed  and  may  vary  with  the  coordinates. 
There  are  two  basic  theoretical  approaches  to  a  solution  of 
the  wave  equation  for  the  specified  boundary  conditions.  The 
first  is  wave  theory,  in  which  the  propagation  is  described  in 
terms  of  characteristic  functions  called  normal  modes,  each  of 


which  is  a  solution  of  the  wave  equation.  The  normal  modes 
are  combined  additively  to  satisfy  the  source  conditions  of 
interest.  The  result  is  a  complicated  mathematical  function 
which,  though  adequate  for  computations  on  a  computer,  by 
itself  gives  little  insight  on  the  distribution  of  the  energy 
of  the  source  in  space.  The  other  basic  approach  to  a 


solution  of  the  wave  equation  is  ray  theory.  The  essence  of 
ray  theory  is  that  the  propagation  path  of  individual  rays  can 
be  traced  by  using  Snells  Law.  The  reader  interested  in  an 
introductory  description  of  wave  and  ray  theory  is  referred  to 


a  book  by  Officer  [Ref.  3].  Snell's  Law,  one  of  the  most 
important  practical  tools  of  ray  theory,  describes  the 
refraction  of  sound  rays  in  a  medium  of  variable  sound  speed. 


As  applied  to  an  ocean  of  horizontal  layers  of  isospeed  water. 
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Snell's  law  states  that  the  grazing  angles  dj,  $2, .  •  •of  a  ray 
at  the  layer  boundaries  are  related  to  the  sound  speeds  Cj, 
C2, . .  •  of  the  layers  by 

cos6i  cosO,  COS63  ^  ^  - 

- -  =  - S  =  - i  =  . . .  =a  constant  for  any  one  ray 

Cj  Cj  C3 

This  expression  enables  a  particular  ray  to  be  traced  out  by 
following  it  through  the  successive  layers  into  which  the 
sound  speed  profile  may  have  been  divided. [Ref .  4] 

A.  RAYMOOE 

The  passive  RAYMODE  model  predicts  the  transmission  loss 
of  sound  through  an  ocean  medium.  It  is  a  Navy  Standard 
Passive  Acoustic,  Range  Independent,  Propagation  Loss  Model. 
The  model  was  developed  at  the  Naval  Underwater  System  Center, 
New  London,  by  Dr.  G.A.  Leibiger  in  1979  [Ref.  5;p.  1]. 

RAYMODE  has  undergone  many  revisions  since  it  was  first 
presented,  but  little  has  been  written  about  the  model  since 
1985.  The  Navy  is  currently  using  Revision  8.1  of  the  RAYMODE 
program  [Ref.  6:p.  iv] .  Although  the  program  has  been  updated 
to  simplify  user  interface  and  to  utilize  new  computer 
technology,  the  basic  physics  of  the  model  for  calculating 
propagation  losses  has  changed  very  little.  The  RAYMODE  model 
combines  the  wave  and  ray  theories,  described  above,  using  a 
normal  mode  summation  technique  for  each  propagation  path. 
RAYMODE  gives  control  of  many  of  the  propagation  loss 
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variables  to  the  user  to  provide  the  most  precise  definition 
of  the  environment.  The  propagation  losses  are  computed  for 
each  path  type  (e.g.,  direct  path,  convergence  zone,  or  bottom 
bounce) ,  at  user-def ined  ranges  and  summed  for  each  source 
angle,  or  for  each  of  a  group  of  normal  modes.  The  total 
field  at  each  range  is  determined  by  summing  the  contributions 
over  all  path  types. 

The  operation  of  the  RAYMODE  model  as  described  by  R.  C. 
Medeiros  in  his  overview  of  the  model  [Ref.  2:pp.  4-12]  is 
summarized  below: 

1.  The  operator  is  required  to  provide  environmental  data 
and  program  control  parameters,  or  indicate  which  historical 
database  to  use.  These  inputs  include: 

a.  Environmental  data: 

Sound  speed  versus  ocean  depth  profile 

Bottom  vi«ipth 

Bottom  type  code  (BLUG) 

Hind  speed 

b.  Geometrical  data: 

Source  depth 
Receiver  depth 

c.  Source  frequency 

d.  Performance  data: 

Minimum  and  maximum  performance  ranges 
Performance  range  step 

2.  Input  data  are  processed  by  the  program  in  preparation 
for  the  actual  computation  and  summation  of  propagation 
losses. 
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a.  Propagation  sound  speeds  are  computed  which 
correspond  to  the  input  source  and  receiver  depths  by  linear 
interpolation  between  the  entered  profile  depths.  These 
depths  and  their  associated  sound  speeds  are  inserted  into  the 
working  profile.  Together,  they  define  the  end-points  of  the 
propagation  paths. 

b.  The  bottom  depth  is  inserted  into  the  working 
profile  and  a  corresponding  sound  speed  is  extrapolated  or 
interpolated.  The  profile  is  then  modified  by  the  program  to 
account  for  the  effect  of  earth  curvature. 

c.  A  gradient  profile  is  computed  from  adjacent  values 
of  c  as  a  function  of  depth.  The  sound  speed  gradient, 
g  =  Lc/Lz,  is  the  ratio  of  the  change  in  sound  speed  and  the 
depth  change  over  which  the  sound  speed  change  occurs. 

d.  The  profile  is  thus  partitioned  into  wave  number 
domains,  (see  Fig.  2)  [Ref.  2:p.  5]  testing  the  profile  for 
localized  maxima  and  minima  which  indicate  those  water  column 
characteristics,  i.e.,  ducts  and  layers,  that  define  the 
propagation  paths.  The  wavenumber 


wavenumber  = 


Frequency ( rad/ sec ) 
soundspeed {yard/ sec) 


(3) 


is  the  parameter  which  most  effectively  specifies  propagation 
path  characteristics  such  as  propagation  angles  and  depths  and 
ranges  associated  with  different  propagation  modes. 
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Where: 

Cl 

Sound  speed  ac 

the 

surface 

CNUP 

■  Sound  speed  ac 

the 

source  depth 

2S 

•  Source  depth 

C2,Z2 

■  Local  ainiaua 

C3.23 

■  Local  Mxiaua 

C4»^4 

■  Local  Baxiausi 

Zb 

•  BoCCoa  depth 

j 

"  Vaveniaiber  interval 

.  counter 

Figure  2  -  Sound  Speed  Profile  Indicating  Criteria  for 
Wavenumber  Partitioning 
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3.  Once  these  wavenumber  domains  are  determined,  they  are 
interpreted  by  the  program  as  angular  intervals,  a  bounded  set 
of  sound  propagation  angles. 

4.  For  each  path,  two  more  quantities  are  defined.  The 
first  is  the  number  of  cycles  that  exist  along  the  path  range 
between  the  source  and  receiver.  A  full  cycle  is  defined  as 
the  path  from  surface  or  apex  to  the  bottom  or  nadir  and 
return.  A  partial  cycle  is  the  path  from  receiver  to  surface 
or  from  the  source  to  surface.  The  total  number  of  cycles  and 
the  partial  cycles  when  appropriately  combined  yield  the  path 
range  for  a  given  propagation  path.  Also  a  phase  change 
occurs  at  the  turning  points  which  must  be  accounted  for.  One 
additional  parameter  to  be  computed  is  the  number  of  modes 
which  exist  for  the  wavenumber  interval  under  examination. 
For  some  paths  there  may  thousands  of  modes,  assuming  the 
frequency  is  sufficiently  high. 

When  the  number  of  modes  propagated  in  an  interval 
becomes  too  large,  the  interval  is  subdivided  into  what  are 
called  K-regions.  A  K-region  is  the  interval  between  two 
wavenumbers  that  defines  a  propagation  path.  Referring  back 
to  Fig.  2,  the  sound  speed  profile  is  partitioned  into 
intervals  bounded  by  a  minimum  and  maximum  sound  speed.  By 
converting  these  sound  speeds  into  wavenumbers  the  interval 
between  the  smallest  and  largest  wavenumber  is  subdivided  into 
a  predetermined  number  of  K-region  intervals.  The  K-region 
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Interval  is  then  converted  into  a  propagation  angle  using 
Snell's  lav. 

Each  K-'region  interval  corresponds  to  an  initial  angle  of 
a  propagated  ray  that  is  extended  down  the  range  path  with 
curvature  based  on  the  application  of  Snell's  law  to  the  sound 
speed  profile.  At  each  range  increment  the  sum  of  all  the 
mode  contributions  from  each  K-region  interval  is  made 
resulting  in  a  single  combined  value  of  acoustic  pressure 
field.  The  individual  contributions  are  summed  at  each  range 
for  all  the  rays  in  an  angular  interval,  which  is  then  summed 
to  those  of  previous  angular  intervals  to  obtain  the  total 
acoustic  pressure  field  at  a  given  range.  After  correction 
for  boundary  losses  and  sonar  beam  deviation  loss  the  result 
is  converted  to  decibels  and  output  as  propagation  loss 
versus  range.  [Ref.  l:pp.  4-6] 

When  the  propagation  angle  increases,  as  it  does  in 
bottom  bounce  regions,  the  number  of  modes  propagated  rises, 
increasing  the  number  of  iterations  necessary.  To  avoid  this, 
when  the  number  of  modes  exceeds  10,  the  model  integrates  over 
a  series  of  K-regions  and  power  sums  or  coherently  sums  the 
contributions  of  each. 

Two  types  of  propagation  losses  are  computed,  coherent 
and  random  incoherent.  Coherent  losses  are  those  attributed 
to  ray  paths  with  an  exact  phase  relationship  between  all 
paths.  Random  losses  assume  a  random  distribution  of  phase. 
Because  coherent  propagation  losses  tend  to  vary  widely,  small 
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modifications  in  input  data  can  cause  large  changes  in  output 
data.  The  operational  RAYMODE  model  utilizes  the  random  phase 
propagation  loss  predictions  outputs  which  tend  to  be 
smoothed  over  the  performance  range  span,  providing  an 
expected  value  of  propagation  loss.  [Ref.  2:p.  12] 

B.  ASTRAL 

ASTRAL  (ASEPS  Transmission  Loss)  is  a  range  dependent, 
low  frequency  propagation  loss  model  developed  by  Science 
Applications  International  Corporation  [Ref.  7:p.ii].  Since 
ASTRAL'S  original  version  in  1979  it  has  undergone  several 
modifications.  The  addition  of  a  surface  duct  and  convergence 
zone  capability  has  increased  its  applicability  and  utility. 
ASTRAL  4.1  is  currently  used  as  one  of  the  Navy  Standard  Range 
Dependent  Propagation  Loss  models  [Ref.  7:p.3].  The  remainder 
of  this  section  is  a  summary  of  the  basic  physics  and  theory 
of  the  ASTRAL  model  as  described  in  Refs.  7  and  8. 

The  calculations  in  ASTRAL  are  carried  out  in  two  parts, 
the  near  field  and  the  far  field.  The  near  field  calculations 
are  based  on  a  ray  trace  scheme.  The  number  of  rays  traced  is 
set  to  50  in  the  operational  version.  The  rays  are  traced  out 
to  a  range  specified  by  the  user.  The  angular  spread  between 
rays  is  adjusted  so  that  80%  of  the  rays  have  launch  angles 
less  than  45  degrees  and  20%  have  launch  angles  greater  than 
45  degrees. 
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In  the  far  field  each  ray  is  associated  with  a  set  of 
modes  or  "smodes**.  These  saodes  approximate  the  envelope  of 
the  oscillatory  depth  function  of  true  normal  modes,  where 
pressure  can  be  expressed  in  the  form 

P  =  ( Zg)  Ua  { Z^)  exp  ( ik^)  ( 4 ) 

where  Z,  and  Z,  are  the  source  and  receiver  depths 
respectively.  The  basic  ASTRAL  model  begins  with  the  normal 
mode  equation  [Ref.  7:p.  12] 

P{r,  2)  *  {z)Uj ZJ  exp  ( jk„r)  < 5 ) 

where  is  a  constant  related  to  the  source  strength  divided 
by  the  square  root  of  the  range  r,  the  U  are  normalized  depth 
functions,  is  a  reference  wave  number  (equal  to  2nf/C0  , 

where  is  a  reference  sound  speed) ,  and  the  asymptotic  form 
has  been  used  for  the  Hankel  function.  The  symbol  i  is  the 
square  root  of  minus  one.  The  intensity  is  assumed  to  be  in 
the  form 

J  =  A^ll„Y^{z^)Y^{z,)/k0  (6) 

where 

Y^{z)  =  [r^(0„)tane„(z)]-i  (7) 

and  where  apex>to-apex  cycle  range  of  a  ray 

launched  with  angle  d„.  The  summation  over  modes  is  replaced 
by  a  summation  over  sets  of  modes,  so  Equation  4  can  be 
%nritten  as 
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(8) 


I  =  (z,)  yj  (z^)  x^. 

where  each  value  of  the  running  index  corresponds  to  a  set  of 
inodes,  or  "smodes."  The  variable  is  given  by  N^/k„,  where 
N.  is  the  number  of  modes  in  the  smode.  In  the  ASTRAL  code, 
these  average  values  are  actually  calculated  at  specific  ray 
angles  measured  relative  to  the  axis  of  minimum  sound  speed 
and  thus  a  smode  corresponds  to  an  angle  bin.  Using  the  phase 
integral  where 

^  flsiniQg^iz)) /c(z)]dz  (9) 

then  .  If  i|r,g/ (sin6d6)  is  assumed  constant,  then  the 

pressure  can  be  expressed  as 

Where 


yjz)  -*  ^yj(z) 


(11) 


and 

X„  =  cos[«oi|r„(z„)]cos[«£,i|f„(z,)]exp{i/v')  .  (12) 

The  smode  is  then  extended  at  the  turning  points  via  Airy 
functions.  The  smodes  are  propagated  in  the  far  field  until 
an  environmental  change  occurs.  If  the  structure  of  the 
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sound-speed  field  changes  then  a  new  smode  assignment  must  be 
made  using  the  new  turning  point  structure  and  assuming 
adiabatic  invariance.  As  the  environment  changes,  the  phase 
integral  is  kept  constant  by  choosing  a  new  turning  point 
sound  speed.  [Ref.  8] 

Convergence  zone  effects  can  be  added  if  each  smode  is 
assumed  to  be  composed  of  four  parts  that  correspond  to  the 
four  different  up/down  arrivals  at  both  source  and  receiver. 
Each  of  the  components  is  propagated  in  range  using  the  cycle 
length  of  the  smode  and  the  intensity  of  each  is  assumed  to  be 
proportional  to  the  square  of  a  sin(x)/x  function.  At  each 
range  the  intensity  is  summed  and  smoothed. [Ref .  8]  The  full 
derivation  of  this  concept  is  beyond  the  scope  of  this  paper. 
For  a  full  description  of  Convergence  Zone  ASTRAL  the  reader 
is  referred  to  Ref.  7,  pp.  14-18. 
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III.  BACKORODIID 


A.  DESCRIPTION  OF  AREAS  USED 

Three  areas  were  chosen  for  the  ASTRAL  and  RAYMODE 
comparisons.  The  areas  were  chosen  based  on  the  new 
regionally  focused  Navy  strategy  stated  in  "From  the  Sea...". 
The  specific  coordinates  for  the  areas  investigated  were 
picked  at  random  and  then  rounded  to  the  nearest  degree  for 
ease  of  data  entry.  Each  of  the  three  areas  are  50  nm  by  50 
nm  boxes. 

The  first  area  is  in  the  Eastern  Mediterranean  centered 
at  33N/033E.  Area  1,  located  between  Egypt  and  Cyprus  as 
shown  in  Fig.  3,  is  typical  of  one  in  which  an  air  ASM 
platform  might  have  to  search  for  a  third  world  submarine. 
The  area  is  close  enough  to  the  coast  to  be  considered 
littoral  but  too  deep  to  really  be  considered  shallow  water. 
The  bottom  contour  and  sound  speed  profile  are  shown  in  Fig  4. 
The  surface  temperature  in  the  area  is  60°  F  ,  there  is  a 
layer  depth  of  350  ft,  the  depth  excess  is  525  fms  and  there 
is  a  sound  channel  at  715  ft. 

Area  2  in  the  Gulf  of  Oman,  is  centered  on  24-30N/059E. 
This  region  has  become  a  fairly  high  interest  area  since  the 
Persian  Gulf  War.  The  specific  area  used  for  the  study  is 
shown  in  Fig.  5.  The  sound  speed  profile  for  the  center  of 
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Figure  3  - 


Location  of  Area  1  in  Eastern 
Mediterranean 


Sound  Speed  Profile  for  Area  1 


Bottom  Profile  of  Area  1 
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Figure  4  -  Sound  Speed  Profile  and  Bottom  Profile  for  Area  1 
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Figure  7  -  Location  of  Area  3  in  South  China  Sea 
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the  area,  and  the  bottom  profile  are  shown  in  Fig.  6.  A 
surface  temperature  of  80°  F  results  in  a  large  negative 
gradient  below  the  layer  depth  of  250  ft.  There  is  no  depth 
excess  (depth  needed  for  convergence  zone  propagation  to 
exist)  in  this  area  and  the  sound  channel  is  fairly  deep  at 
approximately  4000  ft. 

The  third  area  is  in  the  South  China  Sea  at  lON/llOE, 
only  75  miles  off  the  coast  of  Vietnam  (Fig.  7) .  The  sound 
speed  profile  is  similar  to  the  one  in  Area  2  with  a  strong 
negative  gradient  below  the  layer  (Fig.  8)  .  The  surface 
temperature  is  75“  F,  and  the  sound  channel  is  at  2000  ft. 
Similar  to  Area  2,  there  is  no  depth  excess. 

B.  DATA  COLLECTION 

The  ASTRAL,  RAYMODE,  and  PE  models  used  for  this  study 
were  run  inside  a  new  developmental  ASW  Tactical  Decision  Aid 
(ASWTDA  version  2. 1.2.1)  which  has  many  ASW  tactical 
applications.  The  ASWTDA  program  is  run  on  a  SUN  work  station. 
The  databases  used  by  ASWTDA  are  the  Oceanographic  and 
Atmospheric  Master  Library  (OAML)  and  the  GEM  database.  The 
GEM  database  is  a  subset  of  OAML  and  was  created  by  ASWTDA 's 
developing  contractor,  SONALYSTS,  Inc.,  in  order  to  reduce  the 
OAML  database  size  and  the  time  required  to  access  the 
information.  Table  I  lists  the  versions  of  the  models  and 
databases  used  by  ASWTDA.  [Ref.  9] 
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Database  or  Model 

ASWTDA  (V.2.1. 2.1)  I 

LFBL 

2.0 

HFBL 

2.0 

DBDBC 

1.0 

GDEM 

3.0 

GEM 

1.0 

SHIP  NOISE 

4.0 

SEDIMENT 

3.0 

PASSIVE  RAYMODE 

8.1 

PE 

3.0 

CZ  ASTRAL 

2.4 

Table  I  -  Models  end  Databases  Used  by  ASWTDA 


LPBL  and  HFBL  are  low  and  high  frequency  bottom  loss 
databases  I  GOEM  is  an  oceanographic  environmental  database  and 
DBDBC  is  the  bathymetry  database.  The  areas  of  interest  were 
generated  inside  ASWTDA* s  "Environmental**  function  with  20 
sample  points  inside  each  50  nm  by  50  nm  box.  Propagation 
loss  profiles  were  generated  for  the  combinations  of  90  and 
400  ft  receivers,  150  and  450  ft  depth  targets,  and  50  and  400 
Hz  sources.  This  produced  8  range  dependent  and  8  range 
independent  profiles  for  each  area.  The  specific  technique 
and  code  used  to  extract  this  data  is  shown  in  Section  B  of 
Appendix  A.  The  propagation  loss  profiles  were  then  evaluated 
using  a  Matlab  subroutine  to  produce  Median  Detection 
Ranges(MDR)  for  Figures  of  Merit  (FOM*s}  of  65, 70, 75, and  80 
dB.  Three  distributive  sonobuoy  patterns  (4x4,  5x6x5,  2x8) 
were  evaluated  at  each  FOM  and  receiver/ source/ frequency 
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combination.  A  pattern  evaluation  program  provided  by  Wagner 
Associates  called  NESST  was  used  to  produce  probabilities  of 
detection  for  the  three  selected  patterns.  A  description  of 
the  NESST  program  and  how  the  CDP's  were  obtained  is  included 
in  Appendix  B. 

C.  ASSUMPTIONS 

For  the  purpose  of  this  study  several  basic  assumptions 
were  made.  The  first  is  that  the  split-step  PE  propagation 
loss  model  is  sufficiently  accurate  to  serve  as  the  basis  for 
comparisons.  The  target,  for  the  probability  of  detection 
calculations,  is  assumed  to  be  uniformly  distributed  in  a  50 
nm  by  50  nm  box.  In  most  littoral  ASW  situations  the  passive 
detection  ranges  are  relatively  short  (1  to  10  nm) .  Therefore, 
even  though  the  data  for  the  entire  area  was  used  for  analysis 
purposes,  the  profiles  generated  were  only  graphed  out  to  30 
nm.  The  sonobuoys  were  considered  to  be  cookie  cutter 
detectors:  if  the  target  was  inside  the  given  range  of  the 
specific  sonobuoy  it  was  counted  as  a  detection  and  if  it  was 
outside  the  range  it  was  considered  as  undetected. 
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!▼.  MALySZS 


A.  AREA  It  RASTBEH  MBOZTIERAXBAM 

Figures  9’-12  conpare  the  ASTRAL  and  RAYNODE  propagation 
loss  profiles  against  PE  generated  profiles.  Analysis  of  the 
ASTRAL  profiles  shows  such  lower  dB  loss  versus  range  than  in 
the  PE  model  or  the  losses  expected  from  spherical  spreading. 
(The  "x*  marks  on  the  graphs  indicate  spherical  spreading  at 
1,2,  and  3  nm.)  The  RAYMODE  profiles  correlate  very  well  with 
the  spherical  spreading  and  the  PE  profiles,  with  the 
exception  of  the  90ft/150ft/400  Hz  run.  This  lack  of 
correlation  in  the  shallow/ shallow  situation  is  most  likely 
due  to  the  fact  that  both  the  source  and  receiver  are  in  the 
mixed  layer  and  the  400  Hz  source  frequency  is  above  the  150 
Hz  cutoff  frequency  for  a  layer  of  that  depth.  The  RAYNODE 
model  does  not  predict  the  surface  duct  and  therefore  gives 
much  shorter  predicted  ranges. 

Predicted  detection  ranges  were  obtained  from  the 
proploss  profiles  by  extracting  the  range  at  which  the  profile 
crossed  the  given  figure  of  merit.  The  ranges  for  each  FOM 
are  tabulated  in  Appendix  A.  The  average  differences  in  the 
detection  ranges  are  shown  in  Table  II. 

The  calculated  detection  ranges  were  used  to  construct 
the  three  sonobuoy  patterns,  with  a  buoy  spacing  of  1.5  MDR 
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Table  II  -  Area  1  Average  Predicted  Detection  Range 
Differences 
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and  a  row  spacing  of  1.67  MDR  as  prescribed  in  the  ASH 
tactical  manuals.  The  pattern  in  each  scenario  that  generated 
the  largest  ciimulatlve  detection  probability  (CDP)  was  used  in 
Fig.  13  to  compare  the  CDP*s  of  each  model.  (The  first  column 
in  the  table  is  the  run  ID  which  is  described  in  more  detail 
in  Appendix  A.)  It  should  be  noted  that  in  the  case  of 
RAYMODE  vs  PE  (Fig.  14)  at  an  FOM  of  70  the  average  difference 
not  including  the  AISSH  run  is  only  4.0%.  At  FOM's  of  65  and 
70  dB  the  average  difference  in  detection  probability  when 
compared  to  PE  is  fairly  close.  At  an  FOM  of  75  dB  however, 
ASTRAL  vs  PE  has  an  average  difference  of  45.1%  while  RAYMODE 
vs  PE  is  only  7.9%.  The  ranges  at  an  FOM  of  80  dB  were  in 
excess  of  15  nm  and  resulted  in  GDP's  of  100  percent  in  all 
three  models  with  the  exception  of  the  cross  layer  high 
frequency  scenarios  in  the  RAYMODE  model.  (The  80  FOM  results 
were  not  included  in  CDP  comparisons  for  Area  1.) 

No  appreciable  difference  in  computation  time  between 
ASTRAL  and  RAYMODE  was  found  using  a  source  frequency  of 
either  50  or  400  Hz.  Typically  both  models  required 
approximately  3  to  5  seconds  to  generate  proploss  curves  at  50 
Hz.  At  400  Hz  the  ASTRAL  model  took  an  average  of  2  seconds 
longer  to  generate  the  propagation  loss  curve  than  the  RAYMODE 
model  which  generally  took  5  sec. 


400ft/450ft/50  Hz  (c  and  d) 
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Figure  10  - 


Proploss  for  Area  1:  400ft/150ft/400  Hz 
(a  and  b) ,  400ft/150ft/50  Hz  (c  and  d) 
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Figure  12  -  Proploss  for  Area  1:  90ft/150ft/400  Hz 
(a  and  b) ,  90£t/ 150ft/ 50  Hz  (c  and  d) 
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Figure  13  -  Cumulative  Detection  Probability  Comparisons  (%) 
ASTRAL-PE  (left)  and  RAYMODE-PE  (right) 
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B.  MIBB  2t  OULF  OF  OMMI 

The  transmission  losses  in  Area  2  were  much  greater  than 
in  Area  1.  The  propagation  loss  profiles  in  Figs.  14-17  show 
a  greater  than  80dB  loss  in  less  than  five  miles  in  all  but 
the  90ft/150ft/400  Hz  case.  The  rapid  loss  of  sound  energy  in 
such  a  short  range  resulted  in  a  much  higher  correlation 
between  all  three  models.  (The  cyclical  appearance  of  the  50 
Hz  PE  proploss  profiles,  seen  here  in  Area  2  and  even  more 
apparent  in  Area  3 ,  is  examined  further  in  Appendix  B) .  The 
average  differences  in  the  predicted  median  detection  ranges 
are  shown  in  Table  III.  ASTRAL  consistently  predicts  longer 
ranges  than  PE,  while  RAYMODE  generally  predicts  ranges 
shorter  than  PE. 

The  small  differences  in  ranges  similarly  resulted  in  CDP 
differences  (Fig.  18)  usually  less  than  10%.  Just  as  in  Area 
1  computations,  no  significant  difference  in  the  computation 
time  between  ASTRAL  and  RAYMODE  was  noted. 
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Table  III  -  Area  2  Average  Predicted  Detection  Range 
Differences 
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(a  and  b) ,  400ft/ 150ft/50  Hz  (c  and  d) 
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Figure  17  - 


Proploss  for  Area  2:  90ft/150ft/400  Hz 
(a  and  b) ,  90ft/150ft/50  Hz  (c  and  d) 
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Figure  18  -  Cumulative  Detection  Probability  Comparisons  (%) 
ASTRAL-PE  (left)  and  RAYMODE-PE  (right) 
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C.  MtEX  3:  SOUTH  CHINH  SEA 

The  propagation  loss  profiles  in  Area  3  (Figs.  19-22) 
showed  very  large  variations.  In  the  first  two  scenarios 
ASTRAL  and  RAYMODE  predict  ranges  two  tines  greater  than  those 
generated  from  PE  for  a  FOM  of  70  dB.  At  an  FOM  of  75  dB  the 
ranges  are  as  much  as  five  times  greater.  Even  though  both 
models  differ  greatly  from  PE  the  correlation  between  ASTRAL 
and  RAYMODE  is  better  than  in  areas  1  and  2.  In  the  400  Hz 
runs  the  RAYMODE  and  ASTRAL  profiles  have  a  much  higher 
agreement  with  PE  than  in  the  50  Hz  runs.  The  average 
differences  for  just  the  50  Hz  runs  for  ASTRAL  vs  PE  and 
RAYMODE  vs  PE,  at  an  FOM  of  70,  are  2.1  nm  and  2.3  nm 
respectively.  For  the  400  Hz  runs  the  differences  are  only 
1.1  and  .55  nm  respectively.  Similar  relationships  are  true 
for  the  other  FOM's  as  well.  The  reason  for  this  is  not 
exactly  known.  The  overall  differences  in  ranges  are  shown  in 
Table  IV. 

The  differences  between  the  detection  probabilities  (Fig. 
23)  for  all  three  models  appear  to  be  similar  in  this  area. 
Approximately  2%  at  FOM  of  65  dB,  and  13%  at  FOM  of  70  dB. 
The  big  difference  in  this  area  is  that  correlation  between 
ASTRAL  and  RAYMODE  is  better  overall  than  either  ASTRAL  to  PE 
or  RAYMODE  to  PE. 
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Figure  23  -  Cumulative  Detection  Probability  Comparisons  (%) 
ASTRAL-PE  (left)  and  RAYMODE-PE  (right) 
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V.  COMCLUSZOIIS 


The  ASTRAL  propagation  loss  model  consistently  predicted 
much  longer  detection  ranges  than  either  the  P£  or  RAYNODE 
models.  This  could  lead  to  incorrectly  spaced  patterns  and  an 
overly  optimistic  probability  of  detection.  RAYMODE 
generally  predicted  ranges  consistent  with  or  slightly  shorter 
than  PE.  This  would  result  in  a  more  conservative  search  and 
provide  a  lower  chance  of  allowing  a  target  to  slip  through  an 
area  undetected.  At  lower  figures  of  merit  the  difference 
between  models  is  fairly  small.  However,  at  higher  FOM's 
there  are  significant  differences  in  transmission  loss 
predictions. 

No  attempt  was  made  to  compare  the  ASTRAL  or  RAYMODE 
predictions  with  empirically  gathered  data,  thus  no 
conclusions  as  to  which  model  better  reflects  reality  is  made. 

With  the  current  computer  technology  (ASQ  212  computer) 
being  installed  on  ASW  aircraft,  processing  capability  and 
storage  memory  for  new  decision  support  systems  are  limited. 
Considering  this  limitation  combined  with  the  evolution  of 
quieter  submarines  and  shorter  detection  ranges,  installing 
the  more  complex  and  computation  intensive  ASTRAL  program  may 
not  be  the  best  choice.  The  passive  range  independent  RAYMODE 
propagation  loss  model,  as  shown  in  this  report,  will  provide 
adequate  transmission  loss  predictions. 
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JkPPBMDZZ  A 


A.  RUH  ZDSMTZFICATZOM 

Each  propagation  loss  run  was  identified  by  a  consistent 
code  as  follows:  (A  nvim  B  C  D  E},  where  A  represents  Area, 
"num"  is  either  1,  2,  or  3,  B  is  the  receiver  depth  (H*400ft, 
S=90ft) ,  C  is  the  source  depth  (M=450ft,  S«150ft) ,  D  is  the 
source  frequency  (H»400  Hz,  L«50  Hz),  and  E  is  either  RD,  RI, 
or  PE  (RI>=ASTRAL,  RI=RAYMODE,  PE=Parabolic  Equation  Model) . 
For  example  AIMMHRD  is  a  run  in  Area  1  with  receiver  at  400ft, 
source  at  450ft,  frequency  of  400  Hz,  and  was  run  using 
ASTRAL. 

B.  A8WTDA  DATA  COLLBCTZOM 

Data  from  the  ASWTDA  PE  and  ASTRAL  models  was  accessed 
after  the  propagation  loss  plot  was  generated.  The  system 
stored  the  current  graph  data  in  a  file  called  OUTPUT.LOG  in 
the  exec/data/tactical  directory  for  ASTRAL  and  in  a  file 
called  for042.dat  in  the  exec  directory  for  PE.  By  opening  a 
new  window  and  changing  to  the  necessary  directory  the  files 
were  copied  to  a  new  file.  If  this  was  not  done  the  next 
graph  created  would  overwrite  the  previous  information.  The 
ASTRAL  data  was  in  ASCII  format  and  immediately  available. 
The  PE  information  was  in  a  binary  form  and  had  to  be 
conveited  via  a  Fortran  program  listed  in  Appendix  B.  The 
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RAYMODE  data  was  not  stored  for  retrieval  as  in  ASTRAL  or  PE. 
Instead,  using  the  "TDA",  "Range  Independent",  "Proploss", 
"Info"  selections  form  the  menu  bar,  a  tabular  summary  of 
RAYMODE  range  and  propagation  loss  is  displayed.  By  using  the 
mouse  and  left  button  to  highlight  the  vectors,  then  "vi"  to 
open  an  editor  in  an  open  window,  the  center  mouse  button  was 
used  to  copy  the  data  to  a  new  text  file  with  the  appropriate 
name.  [Ref.  7] 

The  proploss  files  were  edited  and  made  into  "m"  files 
for  the  purpose  of  creating  and  analyzing  the  graphs  using  the 
Mat lab  program. 

The  following  tables  contain  the  predicted  Median 
detection  ranges  from  all  three  models  at  each  of  the  four 
Figures  of  Merit  (65,  70,  75,  80  dB) . 
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JUPPENDIX  2 


K.  PB  DATA  C01IVER8Z0H 

The  following  Fortran  program  was  used  to  convert  ASWTDA 

binary  format  (PE  output)  to  ASCII  output.  Load  into  any 

directory  and  run  by  calling  filename. [Ref .  8] 

character*32  infile,  outfile 
character*9  when 

integer*2  ititl(40) ,nr ,  nz,nprof,  nbotm,  nthb,nd 

real  f ,zs,xz(20) ,xr, thmin, thmax,thc,dl,d2 

real  tthb(20) ,rthb(20) ,rplot,tl(20) 

print  *, 'Enter  input  file  (for042.dat  format)  name:' 

read(5, *) infile 

open(unit=l,form='unformatted'  ,status“'old'  ,file=infile) 
print  *, 'Enter  output  file  name:' 
read (5, *) outfile 

open (unita:2 ,  status-' new' ,  f  ile»outf  ile) 
c  READ  THE  HEADER  RECORD  OF  THE  TL  FILE 

read ( 1) ititl, when, f, zs, nr, nz, (xz(i) ,isl,nz) 
c  write  out  four  lines  of  header  information 
print  *, ititl 
%nrite(2,*)  ititl 

print  *,when,f ,zs,nr,nz, (xz(i) ,i=l,nz) 
write(2,*)when,f ,zs,nr,nz, (xz(i) ,i=l,nz) 
print  *, 'Range  TL' 
write ( 2, *) 'Range  TL' 

print  *,' - ' 

write(2,*) ' - ' 

nd=l 

print  *, 'Number  of  depths  is',nd 
c  read  the  range  and  tl  values,  print  them  to  screen  and  c 
to  the  output  file 
do  33  j»l,400 

read(l,end»888,err-999)rplot, (tl(i) , i=l,nd) 
writep,*)rplot,tl(i)  ,i=l,nd) 

33  continue 

888  print  *,'End  of  file  on  .042  file' 
close (1) 
close (2) 

print  *, 'Normal  termination' 
stop 

999  close (1) 


clos«(2) 

print  *,'AbnorBal  temination' 

atop 

and 

B.  MB88T 

The  program  used  to  evaluated  the  emulative  Detection 
Probabilities  of  the  three  patterns  of  interest  was  the  NORDA 
Environmentally  Sensitive  Sonobuoy  Tactics  (NESST)  program, 
provided  by  Wagner  and  Assoc.  NESST  has  two  main 
capabilities.  First,  the  performance  of  an  arbitrary  sonobuoy 
pattern  can  be  evaluated  with  respect  to  a  user  specified 
environment.  Second,  for  a  given  environmental  description  of 
a  region,  NESST  will  develop  a  sonobuoy  pattern  which 
optimally  searches  that  region.  [Ref.  11] 

Only  the  first  of  these  two  capabilities  was  used  for 
this  project.  For  each  of  the  search  regions  a  2x8, 4x4,  and 
5x6x5  sonobuoy  pattern  was  constructed  for  each  FOM,  using  1.5 
times  the  NDR  for  sonobuoy  spacing  and  1.67  times  MDR  for  row 
spacing.  To  evaluate  the  sonobuoy  patterns  NESST  records  the 
location  of  each  buoy  in  the  pattern  in  terms  of  grid 
coordinates.  The  NESST  program  used  the  environment  created 
by  ASWTDA  for  the  50  nm  by  50  nm  regions  as  the  basis  for  the 
grid.  The  patterns  were  centered  on  the  area  of  interest  and 
evaluated  for  probability  of  detection  based  on  the  full  50  x 
50  nm  area.  The  sonobuoys  were  evaluated  as  cookie  cutter 
detectors  with  the  detection  ranges  predicted  by  the  given 
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model  for  the  grid  location.  When  the  full  region  of  interest 
was  100  percent  covered  by  the  cookie  cutter  detectors,  the 
pattern  was  assumed  to  have  a  100  percent  probability  of 
detection. 

C.  CYCLIC  EFFECT  IE  PE  PROFILES 

The  cyclic  appearance  seen  in  the  50  Hz  PE  proploss 
profiles  is  due  to  bottom  reflections.  The  PE  model  generates 
transmission  loss  information  using  a  20"  beamwidth.  In 
areas  where  the  sound  speed  profile  is  strongly  negative,  like 
in  Area  2  and  Area  3,  this  narrow  beamwidth  is  refracted 
downward  and  reflected  off  the  bottom  as  seen  in  Fig.  24. 
When  the  beamwidth  is  increased  to  +/>  40"  (see  Fig.  25)  the 
cyclic  effect  disappears.  Both  the  ASTRAL  and  RAYMODE  models 
use  the  wider  beamwidth  for  their  calculations.  This  seems  to 
indicate  that  in  cases  where  there  is  a  large  amount  of  bottom 
interaction  the  PE  model  is  less  than  ideal.  Since  the 
ranges  examined  in  this  report  are  generally  less  than  10  nm, 
the  information  presented  was  not  affected  by  this  cyclic 
effect. 
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Figure  24  -  Range  Dependent  Ray trace 

with  Beamwidth  of  +/-  20“ 


•  10  «e  to  10 

Range  IPtydt) 


Figure  25  *  Range  Dependent  Raytrace  with 
Beamwidth  of  +/-  40" 
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